The flow of heavy flavor in hydrodynamics 
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The flow of charm is calculated in 2+1 ideal hydrodynamics by introducing the charge of cc pair 
assuming that the number of cc pairs is conserved in relativistic heavy-ion collisions. It is found 
that the mean radial flow velocity of charm quarks is smaller than that of bulk matter by 10~15 
% and the measured V2 of heavy-flavor electrons is reproduced up to = 1.5 GeV/c in Au-|-Au 
coUision at RHIC. The same flow is apphed to regenerated J/tp and its V2 is discussed. 



PACS numbers: 

I. INTRODUCTION 

Heavy flavor is one of the important probes for the 
hot dense nuclear matter created in relativistic heavy- 
ion collisions Unlike the light flavors, heavy flavor 
was expected to be much less thermalized in the matter 
due to the relatively smaller cross section. However, it 
was found that Raa of heavy flavor is much suppressed 
in the intermediate and high transverse momentum re- 
gion and that its V2 is not small 3 ■ This implies that 
charm quarks, which comprise the majority of heavy fla- 
vors produced in heavy-ion collisions, are considerably 
thermalized. On the other hand, charm quarks never 
reach chemical equilibrium due to too small cross sec- 
tion for the annihilation or the creation of the cc pair. In 
fact, the measured number of cc pairs in Au-f-Au collision 
at RHIC is almost proportional to the number of binary 
collisions Q. 

Hydrodynamics, which is based on conservation of 
the energy-momentum and various kind of charges, was 
found to be successful in reproducing the px spectrum 
and the elliptic flow of bulk particles in relativistic heavy- 
ion collisions [H, 0] ■ Additionally, from the fact that the 
number of charm pairs is almost conserved in the flreball 
one can introduce the charge of charm pairs. This 
charge is different from the charge of charm flavor, which 
is zero throughout the space-time of the collision as the 
the same number of anticharm quarks cancels the charm 
charge. Once the charge is introduced, its transport can 
be described in hydrodynamics assuming that the charm 
quark interact strongly with the matter. 

The flow velocity of heavy flavor is a key ingredient in 
the V2 of open charms as well as that of the charmonia. 
In most cases, the flow velocity is assumed same as that 
of the bulk particles or is just parameterized according 
to a given assumption (7H9|. In this study, we attempt 
to calculate the flow velocity of charm quarks in hydro- 
dynamics and use it to obtain the V2 of the open as well 
as hidden charm. 
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The paper is organized as follows: In Sec. |TT1 our 
2-1-1 hydrodynamics simulation is presented and tested 
by comparing the calculated V2 of light hadrons with the 
experimental data at RHIC. We introduce the charge of 
charm pairs into the hydrodynamics in Sec. IIIIl and use 
it to calculate the radial and elliptic flows of charm. In 
Sec. IIVI the elliptic flow of regenerated J/ip is obtained 
with the charm flow. The summary is given in Sec. |V] 
and the details on the semileptonic decay of open charm 
in the Appendix. 

II. 2+1 HYDRODYNAMICS FOR BULK 
PARTICLES 

For our 2+1 hydrodynamics, we use the (r, x, y, rf) co- 
ordinate system defined as follows, 

T = \/W~~z^ , 77 = — In — — . (1) 

'2 t~z ^ ' 

After relativistic heavy-ion collision, hot dense nuclear 
matter is produced between the two receding nuclei. The 
matter is assumed to be thermalized at tq = 0.6 fm/c [lo| . 
The local entropy density at initial time is parameterized 
by 

^ = A|(l-a)^+aneou}, (2) 

where npart(coii) is the number density of participants 
(binary collisions) in the Glauber model, defined by 
"part(coii) = dA^part(coii)/(Torf2:dy); The parameters A and 
a are, respectively, taken as 25.5 and 0.11 to reproduce 
the multiplicity of charged particles at RHIC [III, [12] ■ 

After the initial thermalization, the hot dense nuclear 
matter is assumed to expand according to the hydrody- 
namics. The equations of energy-momentum conserva- 
tion in 2+1 hydrodynamics are given as [Tsl. 

a,(TT"") + a,(TTO-) + a,(rr"^) = -p, 
a^(TT°^) + a^(Tr^^) + a;,(TT^-^) = o, 

dr{TT'^'') + d^{TT^y) + dy{TTyy) = 0, (3) 

where T^^ = (e +p)it^it'^ — p g^^ with m^, e and -p being 
respectively the four- velocity of fiow, energy density and 
pressure of the matter. 
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FIG. 1: the speed of sound squared as a function of energy 
density 



For the equations of state of QGP and hadron gas, 
the quasiparticle model based on lattiee data and the 
resonance gas model are adopted respectively [H, [l6j . 
Whence, a first-order phase transition is assumed and 
the critical temperature Tc is 170 MeV. Fig. [1] shows 
the speed of sound squared, cl = dp/de, as a function of 
energy density, approaches 1 /3 at high energy density, 
which is the value for a free gas, and is zero in the mixed 
phase where pressure is constant in the first-order phase 
transition. 

The equations in Eq. ^ are solved numerically by 
using HLLE algorithm (Godunov-type algorithm) [13- 
[Tgj . The simulation is monitored by entropy conservation 
condition lla. 



dSt, 



drj 



dxdyrsjj 



(4) 



'1/2 



It 



where s is entropy density and 7^ = {1—v'^ —Vy) 
is found that the entropy in our simulation is conserved 
within ± 4% error. 

Fig. [5] (a) shows isothermal lines at the impact param- 
eter 6=9 fm. The temperature for kinetic freeze-out is 
assumed to be at 125 MeV. The lifetimes of the QGP, 
the Mixed and the HG phases are respectively 4.6, 5.4 
and 10.1 fm/c. 

The differential yield of particles of type i in Cooper- 
Frye freeze-out formula is 



dN,, 



1 



dymTdrriTdipp {2n)^ Js 



where (pp is the azimuthal angle of momentum p; E/ 
is the hypersurface at freeze-out temperature; fi{x,p) is 
the phase-space distribution of particles, and d^a^ in 
(r, x, y, rj) coordinate system are, respectively. 



p^^ = I 771T cosh(y - 77), p^, py, —— sinh(y - 77) 




(b) 



FIG. 2: (a) isothermal lines in (j/, r) plane, where y is per- 
pendicular to both the beam direction and the direction of 
impact parameter, and (b) Tf{x,y) at fe = 9 fm in Au-|-Au 
collision at RHIC. 



(6) 



d = I l,-d.j;Tf,-dyTf, \Tf{x,y)dridxdy 



where Tf{x,y) is the invariant time r at freeze-out tem- 
perature on the {x, y) plane, such as shown in Fig. [l](b). 
Substituting Eq. ^ into Eq. © [51], 



dN, 



2g, 



p-d-'aix)Mx,p), (5) XTf 



dymTdniTdipp (27r)'^ 



wltKi 



PT ■ i^Tf Ka 



I nniTjA 



.(7) 



where gi is degeneracy factor. The sign behind summa- 
tion is positive for boson and negative for fermion, i^o(i) 
is modified Bessel function, and V_l = {dx,dy). If Tf does 
not depend on the transverse position Eq. ([7]) is 

equivalent to the blast wave model. 

Two comments are to be made in regard to Eq. ([T]). 
First, if temperature at (t/ -|- e, x, y) is higher than (r/ — 
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e, X, y) with infinitesimal positive s, the square bracket in 
Eq. ([7]) should be multiplied by an additional negative 
sign, because the direction of (fiai, is defined as from 
the higher to lower temperature. Second, if the square 
bracket is negative, it is abandoned. Negative p ■ da^ 
means particle moves into the fireball p^ . 

The V2 of light hadrons from our hydrodynamics are 
compared with experimental data in Fig. [3] Here, V2 of 
particle i is defined as the following: 



V2{pt) 



J d(pp cos{2ipp)dNi/{dymTdmTdipp 
J dippdNi/ {dymTdrriTdipp) 



(8) 



The experimental data include the decay of particles 
like p TT+TT, but ours do not. As in other hydrodynam- 
ics simulations, V2 is reproduced well at low pr including 
the ordering of V2 according to the particle masses. 
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FIG. 3: (color online) V2 of tt^, and proton at b 
fm in hydrodynamics and the experimental data in 20~60 % 
central collision of Au+Au at RHIC 
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III. 2+1 HYDRODYNAMICS FOR HEAVY 
FLAVOR 

The charge conservation in 2+1 hydrodynamics is ex- 
pressed as 



drirn) + d^{Tnv^,) + dy{TnVy) = 0, 



(9) 



where n can be any kind of charge density that is con- 
served. The number of cc pairs produced in relativistic 
heavy-ion collisions, which is proportional to the number 
of binary collisions, is assumed to be conserved during the 
fireball expansion. The justification for this assumption 
is the small cross sections for the creation or the annihila- 
tion of the pair and the overall small number densities for 
these heavy particles Therefore, we can define the 



charge of cc pair. This charge is different from the charge 
of charm fiavor, because the charm net charge is always 
canceled by the same number of the anticharm quarks. 
Now, assuming (anti)charm quarks interact with nuclear 
matter as strongly as light quarks, charm quarks follow 
the bulk flow and Eq. ^ can be used for the transport 
of cc pairs in the nuclear matter. We further assume that 
charm has thermal distribution in each fluid cell, though 
it may fail in high-p^ region. More realistic description 
is possible by solving the Langevin equation [22|, [23| . 

The initial charge density of cc is given by Glauber 
model as following [l6j : 



n(To,a;, y) 



AB 



To 



dz 



dz' pA{s,z)pB{b - s,z'), (10) 



where ct^? is the cross section for cc production in p+p 
collision; A and B are mass numbers of nucleus A and 
B respectively; s and h are the transverse spatial vec- 
tor and the impact parameter vector respectively; Pa(b) 
is the Wood-Saxon distribution function of a nucleon in 
nucleus A(B). The transport of heavy flavor is obtained 
by solving Eq. ([3]) and (|9]) with the initial condition of 
Eq. (Uni). 




FIG. 4: the mean radial distances from the center of fireball 
and the mean radial flow velocities of bulk matter and of 
heavy flavor at 6 = 9 fm in Au+Au collision. 

Fig. S] shows the mean radial distances from the center 
of fireball and the mean radial fiow velocities of the bulk 
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matter and of the heavy flavor at 6 = 9 fm in Au+Au col- 
Hsion. For the bulk matter, the radial distance and the 
radial flow velocity are weighted by entropy density, and 
for heavy flavor, by the charge density of heavy flavor, 
both in the fireball frame. As a result, the mean radial 
distance of heavy flavor is smaller than that of the bulk 
matter by 10'--^15 %. And the radial flow velocity is also 
smaller by similar percentages. The reason is as follow- 
ing: The number of cc pairs is proportional to the number 
of binary collisions while the multiplicity of bulk particles 
mainly to the number of participants, as seen in Eq. ^ 
and (jlOp . The number of binary collisions is more con- 
centrated in the central region of fireball than that of the 
participants. So the mean radial distance of heavy flavor 
from the center is smaller than that of bulk matter. On 
the other hand, the flow velocity is inversely proportional 
to the distance from the center, because pressure gradi- 
ent is high at the border of flreball. Therefore, cc pair is 
less accelerated than bulk particle in hydrodynamics. 

The V2 of open charms is the same as Eq. ([5]) except 
that the differential yield is multiplied by the local charm 
fugacity, 7c(r, a;,j/), which is proportional to the charge 
density of cc in local frame, n{T,x,y)/^_i_{T,x,y): 

J dtpp cos{2ipp) J dxdyjciT,x,y) dN/{dymTdmTdipp) 
J dipp J dxdy^ciT, x,y) dN / (dymTdrriTdtpp) ' 

(11) 

Experimentally the elliptic flow of heavy flavor is mea- 
sured through the electrons produced from the semilep- 
tonic decay of heavy flavor. To reproduce the experi- 
mental data, D mesons are generated according to Eq. 
([7]) with the local charm fugacity being multiplied. It is 
assumed that open charms freeze out at Tc- And then 
61.4 % and 38.6 % of them decay into K + er + and 
K* + e~ + De respectively, assuming that the branching 
ratio of semileptonic decay is 100 %. The energy spectra 
of single electrons in both decay channels are presented 
in Appendix. 

Fig. [S] compares the V2 of heavy-flavor electrons at 6 = 
5.1 fm in hydrodynamics and the experimental data in 
minimum bias [3|. The impact parameter corresponds 
to the mid-centrality weighted by the number of cc pairs 
[2^ 1 . It is seen that the 112 in hydrodynamics reproduces 
the experimental data up to around ^ 1.5 GeV. This 
result implies that the heavy flavor in heavy-ion collision 
at RHIC is considerably thermalizcd at low transverse 
momentum. 



IV. THE ELLIPTIC FLOW OF J/^j 

The V2 of J/V' in relativistic heavy-ion collisions is a 
very interesting quantity revealing how much fraction of 
J/ijj is produced through regeneration Q. Ignoring the 
V2 of primordial J/ ip which seems much smaller than that 
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FIG. 5: the V2 of heavy-flavor electrons at 6 = 5.1 fm in 
hydrodynamics and the experimental data in minimum bias 



of regenerated J/V', the V2 of total J /ij: can be approx- 
imated to be the V2 of regenerated J /ip multiplied by 
its fraction among the total observed yield of J/V'- The 
fraction gives valuable information on the thermal prop- 
erties of J/ '0 as well as of the hot nuclear matter created 
in relativistic heavy-ion collision. Because the J/V' is 
composed of charm and anticharm quark, the differential 
yield of regenerated J in V2 is multiplied by the local 
charm fugacity squared: 



J d(/3p cos(2(/9p) J dxdyj'^{T,x,y) dN/{dymTdmTdipp) 
J dipp J dxdyj'^{T, x,y) dN / (dymTdniTdipp) 

(12) 



In Fig. ini the solid line is the V2 of regenerated J /ip 
at 6 = 9 fm, which corresponds to 20~60 % central col- 
lision. In contrast to our previous study @, where only 
radial component of flow was considered, negative V2 is 
not observed at low pr- The freeze-out temperature for 
J/ip is assumed to be as in the open charm case. The 
flgurc shows that we reproduce the first data point at 
Pt = 1.0 GcV; but overestimate the data at three higher 
Pt- This result could be explained if the measured J/tp 
at these higher pT are mostly primordial ones, which has 
no elliptic flow. In the two-component model, the frac- 
tion of regenerated J/ip is only 20 % at 6 = 9 fm [ll|. 
The dashed line is V2 assuming 20 % of J/ip are regen- 
erated ones and the rest primordial ones. It seems from 
the flgure that the fraction of regenerated J/ip is larger 
than 20 % at pT 1-0 GcV, but smaller than 20 % at 
three higher pT- 



5 




(GeV/c) 

FIG. 6: the 112 of regenerated J/?/j at = 9 fm in hydro- 
dynamics (sohd hne) and 20 % of them (dashed hne). The 
experimental data are taken from [25| . 



V. SUMMARY 

The flow of heavy flavor in heavy-ion collision is stud- 
ied in 2-1-1 hydrodynamics by introducing the charge of 
the cc pair. For this, we assume that the number of heavy 
flavor does not change in fireball but that the heavy flavor 
interacts with the nuclear matter created in relativistic 
heavy-ion collision as strongly as the light flavors. 

It is found that the mean radial distance from the cen- 
ter of fireball and the mean radial flow velocity of heavy 
flavor are less than those of bulk particles by 10^15 %. 
This is so because the creation of heavy flavor is more 
concentrated in the central region than that of bulk par- 
ticles. 

The elliptic flow of heavy-flavor electrons is calculated 
in Monte Carlo method and compared with experimental 
data at RHIC. If heavy flavor particle does not interact 
with the nuclear matter, it has no elliptic flow {v2 = 
0). If the interaction is strong enough, heavy flavor par- 
ticle has the maximum elliptic flow which is given by 
hydrodynamics. It is found that our model reproduces 
experimental data up to around = 1.5 GeV. It im- 
plies that the heavy flavors are considerably thermalized 
at low transverse momentum. 

Finally, the V2 of regenerated J/ip is calculated with 
the same flow of charm. By comparing with experimental 
data, it is found that the fraction of regenerated J/ip 
depends on transverse momentum and the the fraction is 
not small at low px- 



Appendix A: semileptonic decay of open charm 

In this study, two main semileptonic decay channels of 
D meson arc considered: D — K + e~ + D,, and D — > 
K* + e~ + De with the average branching ratios being 
6.09 % and 3.84 % respectively. In the limit of massless 
electron, the differential partial widths are 



dr{D K + e- + Ve) 



dq^d cos 9 

dr{D ^ K* + e- + De) 2/(l-cos6')2 



(Al) 



dq'^d cos 
(1 + cos6')2 



pq 



+ 



\H+{q')\' + sm'9\Ho{q')\'}, (A2) 



where 



{mo + mK'fAi{q^) T 2mDpV{q^) 



1 



m 



D 



,2 2mK'{mD + rriK') 



{ml+ml,)A,iq^)-4p'A2iq') 



In Eq. (jAiP '-^ (|A2p . (7^ and p are, respectively, the four 
momentum squared and the magnitude of three momen- 
tum of W boson in the D meson rest frame. 9 is the 
angle between the three momentum of Icpton and that 
of W boson in the W rest frame, /(g^), Ai{q^), l^(q^) 
and A2{q^) are form factors and are parameterized as the 
foUowings: 



1 - ^iVKolo 



(A3) 



where mpoie = 1.89 GeV for f{q^), and the same value 
is applied to Ai{q^), V[q^) and ^2(9^) for simplicity. 
rv = V{Q)IAi{Q) = 1.62±0.08, and = A2(0)/^i(0) = 
0.83 ± 0.05 taken from the particle data book. 

In D meson rest frame, the momenta of K{K*) and W 
boson are respectively 



p^ = [ \lm]^,K') + P^ 0, 0, - p 



(A4) 



(7p = [viD - ^/toI-^^.j +p2, 0, 0, p y (A5) 

Supposing that the momentum of electron in W boson 
rest frame 
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-, -sm9, -sin6'sin(^, -cos 6*], (A6) 
2 2 2 2 ' 



the energy of the electron in D meson rest frame 

= |7(l-/3cos0) - i(<zo + g.cos0), (A7) 
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FIG. 7: the probability density for the electron in D ^ 
K{K*) + + i>e to have energy Ee in D meson rest frame. 



where /3 = —qz/qa and 7 = qo/q, and 
1 



cos 6* = - 2E^ - mo 

p V 

dcos0= -dE"". 
P 



(A8) 
(A9) 



Eq. (|A1I) with the new variable, ^ is then 
dV{D K + e" + Pe) 



dE^ 



'polo 



'■pole 



1 - 



mo 
2E'' 
mo 



dgV 1/(9') P sin' < 
In ( 1 - 







a 

"^polo 

4 



2E^ 



mo 



m 



polo 



a — m 



polo 



'polo 



(AlO) 



where the upper limit of the integration 



a = 2mniJ 



mo ~ 2E<' 



is given from the condition that sin 6 > 0. The available 
range of E'^ is taken to be from to {mo — mji-)(l + 
rriK /mD) 12. The phase space is the same for D — ;> K* + 
e~ + Vf. after niK is replaced by rriK* ■ 

Fig. [7] shows the probability density for an electron 
from D meson decay to have energy E*^ in the D meson 
rest frame. 
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